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Abstract. Probability distributions for correlation functions of particles interacting via 
random-valued fields are discussed as a novel tool for determining the spectrum of a 
theory. In particular, this method is used to determine the energies of universal A-body 
clusters tied to Efimov trimers, for even N, by investigating the distribution of a correla¬ 
tion function of two particles at unitarity. Using numerical evidence that this distribution 
is log-normal, an analytical prediction for the A-dependence of the A-body binding en¬ 
ergies is made. 

Probability distributions are often studied in the context of Monte Carlo calculations, where a finite 
set of field configurations is used to estimate the value of an observable and its associated statistical 
uncertainty. Previously, it has been argued that the mean and variance of the distribution can often be 
estimated based on physical arguments, in efforts to better understand and control the signal-to-noise 
ratios in statistically noisy calculations in El. Therefore, it seems plausible that this argument may be 
turned on its head, such that knowledge of the distribution can be used to make predictions about the 
theory. As an example, probability distributions of correlation functions in certain cases may prove to 
be useful tools for calculations of the spectrum of the theory. 

The concept of probability distributions arises naturally in lattice calculations. As an example, I 
will consider a Euclidean theory of non-relativistic particles interacting via a two-particle interaction, 
which will be mediated on the lattice by a Hubbard-Stratonovich field, rfi. The discretized action of 
the theory may be written S — Kijj, where, 

= (5r,r'T>p,p' + br,T'-lAp,p.(T) , 

£>P.P' = <5p,p'eP'^2" , Ap,p-(T) = bp,p.+C'/2(p_p,)^p (1) 

Ena. Here t denotes Euclidean time and 0p represents the Eourier transform of the Z 2 -valued aux¬ 
iliary field (p. The limit of unitarity may be seen as a non-trivial UV fixed point in the y6-function 
for the two-particle coupling, corresponding to a fine-tuning of the couplings in the lattice theory. In 
particular, one may define a tower of couplings which depend on powers of the momentum transfer, 
C^(P) - s ' ^inOiniv), up to some cutoff No, and tune these to the unitary point, effectively 
removing the first No orders of the effective range expansion for the scattering phase shift Q. 

In this formulation, the auxiliary fields are chosen to live along the temporal links of the lattice, 
so that, along with the choice of open temporal boundary conditions (valid only for zero temperature 
calculations), the matrix K takes on an upper tridiagonal form. One important consequence of these 
choices is that the determinant, det K is simply the determinant of a product of free kinetic opera¬ 
tors, D, which is independent of rp, and therefore makes no contribution to the probability measure. 
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Figure 1. Histogram of the two-body correlator, C 2 , at large 
Euclidean time, r. Inset: Histogram of the logarithm of the two-body 
correlator at the same r. 
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Furthermore, propagators may be written as iterative operations in time, 

= D-^X{t - - 1 ) ( 2 ) 

with /T '(0) = Z) *, which, for large Euclidean time, act as filters for the ground state, 

V<'F|Z:-'[<^,T]|'P) = y ^ , (3) 

/ i / i r—>oo 

<p n 


where denotes a sum over all held conhgurations <p, and £„ are single particle energies corre¬ 
sponding to the eigenstates \n) that have non-zero overlap onto the wavefunction 'E (the overlap factor 
onto the ground state is denoted by Z)- 

Given a hnite set of held conhgurations, the probability distribution for an observable may be 
sampled by calculating the observable on each conhguration (as is done in Eq.j^prior to performing 
the sum). The values of the observable on each conhguration may be plotted as a histogram, with 
the shape of the histogram representing the distribution. Eor the two-particle correlator tuned to 
unitarity and calculated at large Euclidean time using the lattice theory described above, one hnds the 
histogram shown in Eig. The inset shows the histogram obtained by calculating the logarithm of 
the observable, rather than the observable itself, on each sample. This distribution looks remarkably 
Gaussian, characteristic of an observable which obeys the log-normal distribution. 

While the lattice provides a natural tool for sampling probability distributions, the distributions 
themselves are not simply an artifact of lattice calculations, but may be dehned held-theoretically in 
the continuum. Assuming we have a theory of particles interacting through a field, 0, with probability 
measure e then the probability distribution, P, and characteristic function, <1), for a functional 
Civ[(^], may be defined as 


P{x)oc -x), 0(5) oc I 


/■ 


(4) 


The characteristic function is the Eourier transform of the probability distribution, and is itself a 
partition function for a theory that includes the observable in the Lagrangian, and may be amenable 
to perturbative treatments 0161. In particular, VT(s) = -lnO(s) is the generating functional for the 
cumulants, k, of the distribution. 


W(5) = -2 

«=1 


(/5)« 


(5) 


Thus, if Ca^ = {Cn[(P]) is an A-body correlation function (angle brackets denote an average over the 
field (j) according to the appropriate probability measure for (p) composed of N propagators, then k„ 
may be identihed as a sum of connected diagrams containing nxN external legs. 

In this field theoretic language, the central limit theorem may be described in terms of renormal¬ 
ization group (RG) flow, the idea being that by successively block averaging random numbers drawn 
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from a given distribution characterized by a set of cumulants, k„, one finds that the nth cumulant be¬ 
comes rescaled as k„ —> where Nh is the number of block averaging operations (see, e.g., 

0). Thus, cumulants with n > 2 act as irrelevant operators, vanishing in the limit Nb oo. 

The log-normal distribution follows the same RG arguments, however it applies to products of 
random numbers rather than sums. Thus, the log-normal distribution may be described as an RG 
fixed-point of products of random numbers. As described previously, propagators are given by t 
products of D where the fields (p on each time slice are independent, random fields. Since we 

are interested in the limit t ^ co, some RG argument along these lines may be applicable to this 
case. However, recall that is a matrix, and little is known about products of random matrices 

above two-dimensions 0, so at this time a rigorous mathematical treatment along these lines is not 
possible. 

One may understand the distribution, however, on physical grounds by examining the long-time 
behavior of the moments of the A-body correlation function. In the path integral language, the nth 
moment of the distribution, like the nth cumulant, is given by diagrams containing n x N external 
legs. The moments, however, include both connected and disconnected diagrams, and may be shown 
to correspond to n x A-body correlation functions. As an example, if we begin with a two-body 
correlation function composed, for simplicity of this argument, of single particle propagators >S[0, t], 

C2(t) = <(.S[0,t])^), (6) 

then the variance of the probability distribution of this correlation function is given by the product of 
four propagators, averaged over a set of background fields, 

^ , (7) 

This quantity is, in turn, a four-body correlation function. Following this argument, one may easily 
verify that the Ath moment corresponds to a 2A-body correlation function. 

At this point it is important to note that in order for all moments to correspond to physical cor¬ 
relation functions, the formulation must be such that the particle content specified by the probability 
measure is not at odds with the particle content in the operator. For example, many theories include a 
fermion (or boson) determinant as part of the probability measure, with the power of the determinant 
determining the number of degenerate "flavors" of particle in the theory. As noted previously, deter¬ 
minants in this formulation do not contribute to the probability measure, and the particle content of 
the theory is fully determined by the operator (the set of propagators) for a given moment. In a non- 
relativistic theory this is in a sense natural because at zero temperature and zero chemical potential 
there are no particle/antiparticle pairs arising from the vacuum, so the only information about particle 
content must be encoded in the operators measured. 

Finally, because the long-time behavior of correlation functions is dominated by the ground-state 
energy of the system, the long-time behavior of the Ath moment of the two-particle correlator, Adv)"!"). 
is then 

r-(2N) 

Mn{t) ^ ^ (8) 

r—»oo 

where and Zin are the ground-state energy and wavefunction overlap of the 2A-body system, 
respectively. It should be noted that the set of propagators in all cases are inherently symmetric under 
particle interchange, and will have overlap with both identical bosons and distinguishable particles. 

At unitarity, the energies of A-body bosonic systems have been shown to be tied to the spectrum 
of Efimov trimers nano), whose binding energies are known to obey a discrete scaling symmetry, 
for at least A < 16 ifTTHTSl . One method for showing this is to study the ratio of some A-body 
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Figure 2. Ratio of energies predicted by the log-normal 

distribution (green), compared to direct numerical calculations of the 
energies, (blue) fni . (red) HU), and (gray) 02), as well an ansatz 
(purple band) for the N dependence of the energies based on direct 
numerical calculations I13II14I . The red point shows the only 
calculation involving an excited state, where non-universal effects are 
expected to be smaller than for ground-state calculations. 


binding energy to that of a trimer, while varying the short-distance behavior of the theory. If the ratio 
remains constant, this proves that there is no relevant n-body scale (for 3 < n < N), meaning the 
three-body scale is the only scale determining the physics. This implies that for each trimer there 
must be a corresponding A^-body state with the measured binding energy ratio, due to the discrete 
scale invariance. It also implies that this ratio is universal, in the sense that it is insensitive to any 
UV physics. These ratios have been calculated directly using numerical methods (see data in Fig.[^, 
however, at this time there is still notable disagreement among the various calculations. 

It is now clear from the discussion that the energies of A^-body bound states at unitarity should 
have the form = -unA^, where A* is a three-body scale that may be set using the energy of any 
of the Efimov trimers, and is a universal A-body energy ratio. Thus, the moments of the two-body 
correlation function may be shown to have the following form at large Euclidean time, 

Mn . (9) 

T—>00 

The moments of the log-normal distribution are known to be, 

;V^(^LN) ^ ^ ^(LN) ^ ^ (JQ) 

where and cr are the mean and standard deviation of the distribution of the logarithm of the random 
variable, and on the right I have used the constraint that the first moment is a constant, corresponding to 
a correlation function with zero ground state energy. This form implies that the log-normal distribution 
corresponds to correlation functions whose moments are controlled by a single scale, since correlation 
functions decay as the exponential of an energy scale times Euclidean time. This is precisely what 
is expected from the form of the moments of the two-body correlation function at unitarity, Eq. 
Eurthermore, by comparing the log-normal moments with those of the two-body correlator at unitarity, 
one may predict the A-dependence of the universal constants, a 2 N, giving the following energy ratios; 

( 11 ) 

One should note that the expression above is not expected to hold for arbitrarily large A, because as 
the states become increasingly bound, eventually they will become affected by non-universal cutoff 
effects. Eor the (zero effective range) lattice theory which is used to sample the distribution, the energy 
cutoff imposed by the lattice spacing is found to be a factor of ~ 335 times larger than the energy of 
the lowest Efimov trimer, allowing a large separation of scales in which to find universal behavior. 

The prediction based on the log-normal distribution is plotted in Eig.|^ (lower solid line), along 
with a comparison to other works. The upper solid line is an analytic prediction from MM, using 
data produced from a finite range potential, coupled with a procedure for removing the finite effective 
range effects. This conjecture confirmed the leading A^ dependence predicted by this work, but 
with a significantly different prefactor. The points represent direct numerical calculations of the A- 
body binding energies using potential models with varying dependences on non-universal effects such 
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Figure 3. Cumulant expansion (Eq.|12[l for the two-body correlator at 
unitarity cut off at N^. The green points represent the lattice data on a 
16^ X 1000 lattice, while the blue points are from a mock distribution 
created by expanding around unitarity in such a way as to reproduce 
moments corresponding to the energies found in da. 
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as non-zero effective range and three- and higher-body cutoffs. Note that all of these calculations 
(including this work) are based on ground-state calculations, except for the single red point. For this 
point the binding energy of a 6-body state tied to the first excited Efimov trimer was calculated, and 
is therefore expected to be the least affected by non-universal UV effects. 

Finally, note that for odd numbers of particles the correlator is not positive so an exactly log¬ 
normal distribution would be somewhat unexpected. However, in practice it is found that the distri¬ 
butions for odd-numbered correlation functions are also approximately log-normal with only a tiny 
negative contribution. Because all direct numerical calculations of the N-body binding energies in 
other works show smooth behavior with N, I have interpolated the prediction for even N in this plot. 

To determine quantitatively how close the probability distribution for the two-body correlation 
function at unitarity is to a log-normal distribution, we may turn to the cumulant expansion ina, 
which relates the logarithm of the (field averaged) correlation function to the sum of the cumulants of 
the distribution of In C, 


ln<C) = ^ 

/I=l 


^n(ln C) 
nl 


( 12 ) 


Recall that the cumulants of a distribution act like connected «-point operators with increasing di¬ 
mension. Thus the cumulant expansion may be regarded in an effective field theory context: as a 
distribution is driven toward log-normal, cumulants /('„(ln C) with n > 2 become increasingly irrele¬ 
vant E). For the log-normal distribution, cumulants for n > 2 are exactly zero; any deviation from 
log-normal may be quantified by the relative size of the third and higher cumulants. 

In Fig. 1^ the cumulant expansion for the logarithm of the two-body correlation function at large 
time is shown versus the number of cumulants included in the sum, N^- Comparing the value of the 
cumulant expansion cut off dl - 2 to the asymptotic value, we see that the distribution deviates 
by approximately 2% from the log-normal distribution. This small discrepancy is likely due to sen¬ 
sitivity to lattice artifacts for large moments of C 2 . In the same figure an attempt has been made 
to quantitatively determine how much variations in the energies affect the cumulants of the distribu¬ 
tion. An example of a mock distribution generated by expanding around log-normal and fitting the 
coefficients to produce the energies calculated in lIT^ is shown by the blue points. Using different 
parameterizations of the expansion around log-normal and fitting different combinations of energies, 
the discrepancy from log-normal (given by the value of the cumulant expansion at third order) for 
the energies from 112 varies between ~ 17% - 30%, which is approximately the same as the largest 
discrepancy between the energies from ifT^ and those predicted by an exactly log-normal distribution. 
Using similar expansions around log-normal and allowing for the ~ 2% deviation seen in the lattice 
data results in energies that differ from the log-normal predictions by no more than a few percent. 

I would like to emphasize that the log-normal distribution itself implies that these 2A-body en¬ 
ergies are universal, because all moments depend on a single energy scale, and are therefore not 
particularly sensitive to any non-universal effects. Thus, the lattice data, which shows a log-normal 
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distribution, gives further confidence to the claim that there are universal A^-body states tied to Efimov 
trimers. The log-normal distribution implies other universal properties as well. For example, if the 
two-body correlation function is log-normal, then by extension the correlation functions for all 2N- 
body Efimov states must also be log-normal, but with rescaled parameters fi and cr^. The logarithm of 


the nth moment of a 2A^-body correlator corresponds to = iNn{Nn - l)E''yT. Comparing this 


7(4) n 


to the moments of the log-normal distribution (Eq. 10 1 , gives /r = \NE'^'T ,cr^ — 


-N^E^^h for the 

distribution of the 2A^-body correlation function. Thus, the distribution itself may be seen as universal 
for all correlation functions of states tied to Efimov trimers. 

It is compelling that as the beta function for the two-particle coupling approaches an RG fixed 
point, the distribution for the two-particle correlator appears to as well. It may be interesting to 
explore an effective field theory type (or possibly some perturbative) expansion around the log-normal 
distribution mi, much in the way that the KS W expansion lfT9l IMI for nuclear effective field theory 
expands around the unitary fixed point, either analytically using Eq. [^or numerically using lattice 
data. Such a program, in addition to giving further insight into the connection between the log-normal 
distribution and unitarity, would allow calculations of non-universal effects at large N, as well as any 
possible deviations from log-normal within the universal regime. 
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